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Using fourth-order perturbation theory, a general formula for the van der Waals potential of two 
neutral, unpolarized, ground-state atoms in the presence of an arbitrary arrangement of dispersing 
and absorbing magnetodielectric bodies is derived. The theory is applied to two atoms in bulk 
material and in front of a planar multilayer system, with special emphasis on the cases of a perfectly 
reflecting plate and a semi-infinite half space. It is demonstrated that the enhancement and reduction 
of the two-atom interaction due to the presence of a perfectly reflecting plate can be understood, 
at least in the nonretarded limit, by using the method of image charges. For the semi-infinite half 
space, both analytical and numerical results are presented. 

PACS numbers: 12.20.-m, 42.50.Vk, 34.20.-b, 42.50.Nn 



I. INTRODUCTION 

The dispersive interaction between two neutral, un- 
polarized, ground-state atoms — commonly known as the 
van der Waals (vdW) interaction — may be regarded, in 
the nonretarded limit, i.e., for small interatomic separa- 
tions, as the mutual interaction of the fluctuating electric 
dipole moments of the atoms in the ground state. It was 
first calculated in this limit by London using perturba- 
tion theory, the leading-order result being an attractive 
potential proportional to r -6 , where r denotes the in- 
teratomic separation [1]. In the retarded limit, i.e., for 
large interatomic separations, the interaction is due to 
the ground-state fluctuations of both the atomic dipole 
moments and the electromagnetic far field. This was first 
demonstrated by Casimir and Polder, who identified the 
vdW interaction as the position-dependent shift of the 
system's ground-state energy due to the coupling be- 
tween the atoms and the electromagnetic field [2] . Using 
a normal-mode expansion of the electromagnetic field and 
calculating the energy shift in leading-order perturba- 
tion theory, they generalized the (nonretarded) London 
potential to arbitrary distances between the two atoms, 
where in particular in the retarded limit the potential 
was found to vary as r~ 7 . 

The theory has been extended in many respects, and 
various factors affecting the vdW interaction have been 
studied. Based on a calculation of photon scattering am- 
plitudes, Feinberg and Sucher extended the theory to 
magnetically polarizable atoms [3]. They found that the 
vdW interaction of two magnetically polarizable atoms is 
again attractive, while for two atoms of opposed type — 
one magnetically and one electrically polarizable — a re- 
pulsive vdW force may be observed. Later on, it was 
demonstrated that in the case of two atoms of opposed 
type the nonretarded potential is proportional to r , in 
contrast to the r _6 -dependence of the nonretarded poten- 
tial of equal-type atoms [4]. The Feinberg-Sucher result 



was extended to particles exhibiting crossed polarizabil- 
ities [5]. Further studies have also included the cases of 
one [6] or both atoms [7, 8] being excited, leading to po- 
tentials that vary as r~ 6 and r~ 2 in the nonretarded and 
retarded limits, respectively. Thermal photons present 
for any nonzero temperature have been shown to lead, in 
the retarded limit, to a change of the vdW potential of 
two ground-state atoms from a r~~' ' - to a r~ 6 -dependencc 
as soon as the interatomic separation exceeds the wave- 
length of the dominant photons [9-12]. Modifications 
of the vdW interaction due to external fields have been 
shown to lead to a potential varying as r -3 in the nonre- 
tarded limit when the applied field is unidirectional [13]. 
Generalizations of the vdW interaction to the three- [14- 
17] and JV-atom case [18, 19] were addressed first in the 
nonretarded limit and later for arbitrary interatomic sep- 
arations, where the potentials were seen to depend on the 
relative positions of the atoms in a rather complicated 
way. 

Van der Waals interactions play an important role in 
the understanding of many phenomena — mostly in the 
field of surface science, such as surface tension [20, 21], 
adhesion [22], and capillarity [23], but also in chemical 
physics, such as colloidal interactions [20, 24] and sta- 
bility [25]. However, application of the theoretical re- 
sults to these phenomena requires taking into account 
the influence of media on the atom-atom interaction. An 
expression for the vdW interaction of two ground-state 
atoms in the presence of dielectric media was first ob- 
tained by Mahanty and Ninham based on a semiclas- 
sical approach [26-28], and was applied to the case of 
two atoms placed between two planar, perfectly conduct- 
ing plates [27] . The situation of two atoms between two 
perfectly conducting plates was later reconsidered taking 
into account finite temperature effects [29] . Other scenar- 
ios such as two atoms placed within a planar dielectric 
three-layer geometry [30] or two anisotropic molecules in 
front of a dielectric half space or within a planar dielectric 
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cavity have also been studied [31]. 

In this paper we present an exact derivation of a very 
general formula for the vdW potential of two ground- 
state atoms in the presence of an arbitrary arrangement 
of dispersing and absorbing magnetodielectric bodies. 
Based on macroscopic quantum electrodynamics in lin- 
early, locally and causally responding media, and starting 
from the multipolar coupling Hamiltonian for the atom- 
field interaction in electric-dipolc approximation, we cal- 
culate the vdW potential in leading, fourth-order per- 
turbation theory. We then apply the general result to 
the cases that the two atoms are placed (i) within bulk 
material and (ii) in front of a planar magnetodielectric 
multilayer system. 

The paper is organized as follows. In Sec. II the atom- 
field interaction Hamiltonian in its multipolar coupling 
form is presented. The derivation of the general formula 
for the vdW potential is given in Sec. Ill, and Sec. IV 
is devoted to the applications mentioned, where a de- 
tailed analytical as well as numerical analysis is given. 
Finally, the paper ends with a summary and conclusions 
in Sec. V. 



II. MULTIPOLAR-COUPLING HAMILTONIAN 

The Hamiltonian for a system consisting of nonrcla- 
tivistic charged particles a (each particle having charge 
q a , mass m a , position r a , and canonically conjugate mo- 
mentum p a ) interacting with the electromagnetic field in 
the presence of dispersing and absorbing magnetodielec- 
tric bodies is given by [32, 33] 

^= E /d 3 r^°°d W ^fl(r,^).f A (r, W ) 



X—e,m 

+E 



2m G 



p a - g Q A(f a ) 



+ 



+ 



where 



and 



J d 3 r/5 p (r)^(r), 

a 

Pp(r') 



<p p (r) = / dV 



47T£o r — r' 



d 3 rpp(r)<£p(r) 
(1) 

(2) 
(3) 



are the charge density and scalar potential of the parti- 
cles, respectively. The Bosonic fields f A (r, uj) and fj(r, uj) 
are the canonically conjugate variables that describe the 
combined system of the electromagnetic field and the (in- 
homogeneous) magnetodielectric medium, including the 
dissipative system responsible for absorption, 

; Ai (r,w),/l v (r>')l =5 xx ,5 w 5{r-v')5{u-J), (4) 



(5) 



where A = e (X — m) refers to the electric (magnetic) exci- 
tations. The vector potential A(r) and the scalar poten- 
tial ip(r) of the medium-assisted electromagnetic field can 
in Coulomb gauge be expressed in terms of the dynamical 
variables f\(r,u>) and f^(r, u>) as 



l« = r 

Jo 



V£(r) 



dw^w^ET^u;) +H.c, 



dwEH(r,w) +H.c, 



with 



E(r,w) 
where 



]T / d 3 r'G A (r,r».f A (r', W ), 

A— e,m * 



G e (r,r',u) = i 
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TT£q 



Im£(r', uj) G(r, r', uj), 



(6) 
(7) 

(8) 
(9) 



G m (r, r', uj) = — i — G(r, r', uj) xV-m/ — — Imntr 1 , uj), 
c V tt£o 

(10) 

[G(r,r',u))x'V r >~\. : .=ejkid' l Gik(r,r',u)), and _L (||) denotes 
transverse (longitudinal) vector fields. In Eqs. (9) and 
(10), G(r, r', uj) is the classical Green tensor obeying the 
equation 



VxK(r,w)Vx- 



■e(r,ui) 



G(r,r» = 5(r-r') (11) 



together with the boundary condition at infinity. All rel- 
evant characteristics of the macroscopic bodies enter the 
theory via the space- and frequency-dependent complex 
permittivity e{r,u>) and permeability fi(r,uj) = K~ 1 (r,uj), 
with the real and imaginary parts of e(r, uj) and n(r, uj) 
satisfying the Kramers-Kronig relations. Note that the 
Green tensor obeys the useful properties [32] 



G*(r,r',w) = G(ry,-u*), 
G(r,r» = G T (r',r,^), 

J2 I d 3 S G A (r,s,^-G+(r',s,c) 
\ J 



A— e,m ' 



^A*0 2t /-f/ I \ 

uj ImG(r,r , w). 

7T 



(12) 
(13) 



(14) 



If the charged particles constitute a system of neutral 
atoms and/or molecules (briefly referred to as atoms in 
the following) labelled by A, ^2 a eA 1 a = ^ tncn ^ ^ s 
convenient to employ the Hamiltonian in the multipolar- 
coupling form, which can be obtained from the minimal- 
coupling form (1) via a Power-Zienau transformation 



U = exp 



- f d 3 r^P A (r)-A(r) 



(15) 
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where the polarization of atom A is given by 

Pa(t) = ^2 1ar a / d\S(r - r A - Ar Q ), (16) 



with 



(17) 



denoting the particle coordinates relative to the center of 
mass 



r A = > 

^ m A 



(18) 



of atom A [m A = J2aeA m a)- We assume that all the 
atoms are (i) essentially at rest, m a /mA — > 0, (ii) small 
compared to the wavelength of the relevant field compo- 
nents, r a ^>r A , and (iii) well separated from each other, 

J d 3 rP A (r) • P B (r) = S AB J d 3 rP 2 A (r). (19) 

Under these assumptions, the Hamiltonian in the mul- 
tipolar coupling scheme can be obtained from Eqs. (1) 
and (15) in complete analogy to the procedure outlined 
in Ref. [34], resulting in 



H = H F + H A + ^2 Hafi 



(20) 



where 



H F = V d 3 r du;^f A t(r, W )-f A (r,w), (21) 

X=e,m J J ° 



aeA 



H AF = -d A ■ E(r A )+J2 y^t> a ■ [r a * B(r A )} 



a£A 



2m c 



a£A 



In Eq. (23), 



aeA aeA 



(23) 



(24) 



is the electric dipole moment of atom A, and the electric 
and induction fields are given by 

/•OO 

E(r) = / dwE(r,w) + H.c, (25) 
Jo 

with E(r, u) from Eq. (8), and 

/•OO 

B(r) = / dwB(r,w)+H.c, (26) 
Jo 

B(r, lu) = (iuj)- 1 V x E(r, u). (27) 



Note that in the multipolar-coupling scheme E(r) has the 
physical meaning of a displacement field w.r.t. the polar- 
ization of the atoms. Finally, in the case of atoms which 
are not magnetically polarizable, we may omit the sec- 
ond and third terms in Eq. (23) so that Eq. (23) reduces 
to the well-known electric-dipole term 

Haf = -d A • E(f A ). (28) 



III. THE VAN DER WAALS POTENTIAL 

Let us consider two neutral, ground-state atoms A 
and B at given positions r A and rg in the presence of 
arbitrarily shaped magnctodiclcctric bodies. Denoting 
by \n A (B)) the (unperturbed) energy eigenstates of atom 
A(B) 1 we may represent the atomic Hamiltonian H A ( B ), 
Eq. (22), in the form 



H A (B) - ^2 E A{B)\ n A(B)){n A ( B )\. 



(29) 



Restricting our attention to the electric-dipole approxi- 
mation, the interaction Hamiltonian H A ( B )f reads, ac- 
cording to Eq. (28) [r A{B ) ^ ^A(B)], 



H 



A{B)F 



-EEWKb)!^) -Mr a(b)), (30) 



where d A ™ B) = {n A (B)\&A(B)\rn A ( B )), and E(r) is given by 
Eq. (25) together with Eq. (8). Further, let |{0}), |l( Q )), 
and |1^, l*- 7 ') be the vacuum, single-, and two-quantum 
excited states of the combined system consisting of the 
electromagnetic field and the bodies, respectively, 



Ai(r,w)|{0}> = 0, 
/l <a (r a , WQ )|{0}) = |lW), 



(31) 
(32) 



^/U^,^)/^^,^)!^})^!!^,!^) (33) 

[the corresponding single- und two-excitation energies are 
respectively hw a and h(cup + w 7 )]. 

Following Casimir's and Polder's approach [2] (see also 
Ref. [35]), we identify the two- atom vdW interaction with 
the position-dependent shift of the ground-state energy 
AE A b calculated in leading-order perturbation theory 
according to 



AE 



ab 



E 



\Haf+H bf \III}{III\Haf + Hbf\II} 



i, ii, hi 



(Ej — Eq) 

{II\H af +Hbf\I)(I\H af + H B f\0) 
(En — E )(Ejij — Eq) 



, (34) 



where the primed sum indicates that only intermediate 
states | J), | II), and \III) other than the (unperturbed) 
ground state of the overall system, 



|0) = |<U>|0 B >|{0}>, 



(35) 
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are included in the summations. Note that the summa- 
tions include position and frequency integrals. 

From Eq. (30), by considering only two-atom virtual 
processes, it can be inferred that the intermediate states 
|7) and \III) have one of the atoms excited and one body- 
assisted field excitation present, while the intermediate 
states 1 77) can be of three types: (i) both atoms in the 
ground state with two field excitations present, (ii) both 



atoms excited with no field excitation present, and (iii) 
both atoms excited with two field excitations present. All 
possible intermediate states together with the respective 
energy denominators are listed in Tab. II in App. A. 

Let us consider, e.g., case (1) in this table. Substi- 
tuting the corresponding matrix elements (Al)-(A4) as 
given in App. A into Eq. (34), we derive the contribution 
A£xb{i) to the two-atom energy shift AE A b to be 



AE 



AB(1) 



EE E 



2h 3 



n.ni ii,22,*3,«4 Ai,A2,A3,A4 



n 

j=i 



d 3 r 



didj 



D nm (u)i,U)2, Uls, UJ4) 



+ [df-G\ l {v Al v 1 ^ 1 )] ii [d°" • G\ 3 (r A , r 3 , w 3 )] l3 [dg° ■ G\ 2 (r B , r 2 , c 2 )] ^ [d° B m ■ G\ 4 (r B , r 4 , W4 )] u 6 (3 ^ 
+ [< ■ G* Xl (r A , ri , Wl )] h [d° A n ■ G\ 2 (r A , r 2 , a*)] h [d£° ■ G\ 2 (r B , r 2 , u 2 )] ^ [d° B m ■ G\ 4 (r B , r 4 , W4 )] u 5^ 5 (34) 
+ ^0^,^,^]^^^ (36) 

I 



where 



and 



<5 (Q,3) = ^i^A a A^(r Q - r p )S(u a - 143) (37) 



79„ m (wi,cj 2 , w 3 , uj 4 ) = (w^+w 1 )(w 2 +w 3 )(w5+w 4 ) (38) 



J A(1 



(E n A 



Ml 



simplify Eq. (36) to 



- E A ( B j)/h]. Recalling Eq. (14), we may 



AE 



AB(l) 

1 



Mo \ ~> 

fiTT 2 



dw'wV 2 



X 



D +— j[d A n -ImG(r A ,r B ,uj)-d B m } 

x [d° A n -lmG(r A ,r B ,u; r )-d B m ], (39) 

where 7); and Da are respectively the first and the second 
denominators in Tab. II, and without loss of generality 
we have assumed that the matrix elements of the electric- 
dipole operators are real. 

The contributions AE AB ^ to AE AB which corre- 
spond to the cases (2)-(10) in Table II in App. A can 
be calculated analogously. It turns out that they differ 
from Eq. (39) only in the energy denominators. It is 
not difficult to prove that summation of the energy de- 
nominators under the double frequency integral leads to 
(App. B) 



Xll 



1 



4(u/| +u™ + lu) 



1 

id + u' u — u' 



1 



(40) 



Hence, the two-atom contributions AE AB ^ to the 
fourth-order energy shift lead to the vdW potential 
U AB (r A , r B ) = J2k=i ^ E AB(k) as follows: 

u AB (v A ,vb)=-^y.^^k 

^ 2 tt U n A + Jo Jo 



n.m 



J A I" ^B JO 
:2 ,/2 1 .11 j_ , ,m j_ / ^ 1 



id 2 u' I {id\ + id B 1 + LO) 
(wj + w)K + id) V^ + ^' oj-oj' 

x [d° A n • ImG(r, , r B , w) • d° B m ] [d° A n ■ ImG(r A , r B , J) • dg»] . 

(41) 

To perform the integral over u)', we first use the identity 
ImG= (G- G*)/(2i) and the relation (12) to write 



do/ . 

\ld + Id Id — W 

~ 2~i 



dto' 



7 j ^ImG(r A ,r B , W ') 
! 1 



U> + Id Id — Id' 



7 )uj' 2 G(r A ,r B ,Lo'), 

(42) 



where the poles at u' = —uj and u' = id are to be treated 
as principal values. The Green tensor is analytic in the 
upper half of the complex frequency plane including the 
real axis, apart from a possible pole at the origin. In 
addition, w' 2 G(r A , r B , oj') is well-behaved for vanishing 
uj' [32] . We may therefore replace the integral on the right 
hand side of Eq. (42) by contour integrals along infinitely 
small half-circles surrounding ±u>, and an infinitely large 
half-circle in the upper complex half-plane. The integral 
along the infinitely large half-circle vanishes because [32] 



M 



lim u> 2 G(r A ,r B ,u>) 



0. 



(43) 
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Collecting the contributions from the infinitely small where we have again made use of the relation (12). Sub- 
half-circles, we end up with stitution of Eq. (44) into Eq. (41) leads to 



f 

Jo 



/ i i \ 
da/ —— u' 2 lmG{v Al r B ,uj') 

\U> + OJ U) — U)' J 

= ^L0 2 [G(r A ,v B ,uj) + G*(r A ,r B ,uj)], (44) 



a 2 i r°° 



{[d A n -G(v A ,r B ,w)-d B m ] 2 [d° A n -G*(r A ,r B ,u;)-d B m } 2 } 



= -M. V' — 



1 



K+c)(^+c) + 



[d^.G(r A ,r B)W )-d° B m ] 2 . (45) 



This equation can be further simplified by again using 
contour-integral techniques. It can be seen that the in- 
tegrand in the first integral in Eq. (45) is analytic in the 
first quadrant of the complex frequency plane, including 
the positive real axis. Therefore, it can be replaced by 
contour integrals along an infinitely large quarter-circle 
in the first quadrant and along the positive imaginary 
axis, introducing a purely imaginary frequency, lu = iu. 
The integral along the infinitely large quarter-circle van- 
ishes because of Eq. (43). In a similar way, the second 
integral in Eq. (45) can also be transformed to one over 
the imaginary axis. Combining the contributions from 
the two integrals leads to 



du U^LU^LOg 



[(^ A ) 2 + u 2 ][(^) 2 + u 2 } 



x [df-G(r A ,r B ,iu)-d!g- 



0m]2 



(46) 



An expression of this type was first given in Ref. [28] on 
the basis of a heuristic generalization of the respective 
free-space result. 

Noting that the (lowest-order) atomic ground-state po- 
larizability tensor is (see, e.g., [36]) 



a A(B) ( W )= lim 



, ,n jOii r]nO 

A(B) A(B) A(B) 

r,^o+ h ^ {u A(B) ) 2 -OJ 2 - irjuj 



(47) 



we may rewrite Eq. (46) as 



U AB (v A ,v B ) = [°°du 
27T Jo 



x Tr[ot A (iu)-G(r A ,r B ,iu) ■at B (iu)-G(r B ,r A ,iu)], 

(48) 

where we have used Eq. (13). In particular for atoms, 



which are spherically symmetric, 
ctA(B)(w) = a A(B) {uj)I 



= lim — 

r f ^0+ 3ft 



, ,n \Avn 12 r 
W A(B)l a A(B)l J 



Eq. (48) becomes 

U AB (r A ,r B ) = [ duu 4 a A (iu)a B (iu) 

27T Jo 

x Tr[G(r A ,r B ,iu)-G(r B ,r A ,iu)] 



(49) 



(50) 



The total force acting on atom A and B can be derived 
from the potential 

U(y a ,t b ) = U a (t a ) + U b (y b ) + U ab (v a ,t b ) (51) 

according to 



F A (B) = -V rAiB ^U(r A ,r B ), 



(52) 



where Ua(b) is the single-atom potential (see, e.g., 
Ref. [34]) 

Ua(B){*a(b)) = -^T J Auy2 

x Tr [a A(B) (iu) • G (1) (r A(B) , r A(B) , iu)] , (53) 

with G^ being the scattering part of the Green tensor, 

G(r,r',iu) = G (0) (r,r',m) + G (1) (r,r',zu) (54) 

[G(°\ bulk part]. In particular, the body-assisted force 
acting on atom A(B) due to the presence of atom B(A) 
reads 



' AB(BA) 



U AB (r A ,r B ). 



(55) 



Note that F AB ^ —F BA in general, due to the presence 
of the bodies. 
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IV. APPLICATIONS 
A. Bulk material 

Let us first consider the simplest configuration where 
the two atoms are embedded in a bulk magnetodielectric 
material whose Green tensor reads [33] 



G(r,r',m) = G (0) (r,r',m) 
n{iu) 



Att\y - 



(r-r')(r- r') 



r — r ( 



l\2 



x e 



-n{iu)\r— r \u/c 



where n(iu) = ^s(iu)fi(iu) and 



f(x) = 1 + x + x 2 , 
g(x) = 1 + 3x + 3x 2 , 
£ = c[n(iu)\r — v'\u\~ ■ 



(56) 

(57) 
(58) 
(59) 



Combining Eq. (51) [together with Eqs. (50) and (53)] 
with Eq. (56), we find that (I = \r A — i\b|) 

U(r A ,r B ) = U AB (r A ,r B ) 



h/4 



s: 



du 



^faWM -2n(iu)ul/c 



16tt 3 Z 6 Jo ~" e 2 (iu) 

x 1 3 + 6n(iu)ul/c + 5[n(iu)ul / c] 2 

+ 2[n(iu)ul/c} 3 + [n(iu)ul/c] 4 Y (60) 

which generalizes earlier results [2] on the two-atom vdW 
interaction in free space. Note that in Eq. (60) local-field 
corrections are disregarded. They could be taken into 
account in a similar way as in the case of single-atom 
systems (see, e.g., Ref. [33, 37, 38]). 

In the retarded limit, where I 3> c/w m i n [w m ; n = 
min({o>^,, lo v \A! — A,B; n, v = 1, 2, . . .}), with ld v denot- 
ing the resonance frequencies of the medium] , due to the 
presence of the exponential in the integrand in Eq. (60), 
only small values of u significantly contribute. Hence 
we may approximately replace the atomic polarizabilities 
and the permittivity and permeability of the medium by 
their respective static values, 

aA(B)(iu) - a A (B)(0), e(iu)~e(0), fi(iu) ~ //(0), 

(61) 

and perform the integral in closed form to yield 



U(r A ,r B ) = - 



where 



23hc aA(0)a B (0) 



64tt 3 £2 „(0) e 2( ) 



(62) 



(63) 



Equation (62) reveals that the potential behaves like l~ 7 
just as in the free-space case, but with the coefficient 
being reduced by a factor of [n(0)£ 2 (0)] _1 . 



In the nonretarded limit, where I <C c/[n(0)w max ] [w max 
= max({w^, , u> v \A' = A,B;n,v — 1,2,...})], the inte- 
gral in Eq. (60) is effectively limited to a region where 
e -2n(iu)ui/c ^ j anc j term in curly brackets is approx- 
imately equal to 3, so that 



U(r A ,r B ) = 



where 



Cnr — 



2>h 



lQ7T 3 e 2 



du 



a A (iu)a B (iu) 
e 2 {iu) 



(64) 



(65) 



which shows the Z -6 -dependence also known from the 
free-space case. According to Eq. (60) and Eqs. (62)- 
(65), a bulk magnetodielectric medium tends to inhibit 
the interaction between the atoms, thereby reducing the 
interatomic dispersion force. 



B. Multilayer systems 

Now let the two atoms be in front of a planar mag- 
netodielectric multilayer system consisting of N adjoined 
layers labeled by j (j = 0, 1, 2, .., N— 1) with thicknesses dj 
(rfo^oo), permittivities Sj(uj), and permeabilities (j,j(u>), 
as sketched in Fig. 1. The z axis is perpendicular to 
the layers, with the origin being on the interface between 
layer j = N—l and the free-space region, which can be re- 
garded as layer j = N (d^ — > oo, Sn(^) = 1, = !)■ 
With the coordinate system chosen such that the two 
atoms (in the free-space region) lie in the xz plane, the 
nonzero elements of the scattering part (r A , r B , iu) 
of the Green tensor G(r A ,r B ,iu) in Eq. (50) can be given 



by (App. C) 



1 f°° 

G { £ {yy) (r A ,r B ,iu) = — ^ dqqe 



bNZ~ 



MqX)P ) J 2 (qX) s b N [J (qX)^ ) J 2 (qX)} p 
r N ~ r 



b N 



k 2 



N 



(66) 



G^ zx) (r A ,r B ,iu) = fti- [dqq 2 e-^ J -l^lrl, 

(67) 



N ' 



(68) 



where Z+ = z A + z B , X = x B — x A , J v {x) denotes Bessel 
functions, and 



bj = bj{q, u) = \ —Sj(iu)nj(iu) + q 2 , 



(69) 



kj = kj(q,u) = \Jej(iu)fij(iu) — = rij(iu) — . (70) 
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FIG. 1: Sketch of the planar multilayer medium. 



The (generalized) reflection cocfiicicnts rj with respect 
to the left boundary of the jth layer (j = 1, 2, 3, . . . , N) 
can be obtained from the recurrence relation 



rj = rj(«,u) 



-2fcj_idj_i„(T 



J'-l 



-26j_i i ^tr 



(71) 



j'-i 



'o — 



= [a = s,p), where A| and A^ 1 stand for nj(in) and 
Ej(iu), respectively. 

According to the decomposition (54) of the Green ten- 
sor, the two-atom potential U A b, Eq. (50), can be de- 
composed into three parts, 



U ab(? Ai r s) = U A °l(r A ,r B ) + U ( ab (y Ai y b ) 



(i) 



+U A 2) B (r A ,r B ), 



(72) 



where 



U Ab( v A,Yb) = ( duu 4 a A (iu)a B (iu) 
2n J a 

x Tr [G<°> (r^ , r B , in) • G^ (r B , r A , in)] (73) 

is the bulk-part contribution, which is given by Eq. (60) 
with n(iu) = 1 = n(iu), 

U A l(r A ,r B ) = --^ f duu 4 a A (iu)a B (iu) 
?r Jo 

x Tr [G?(°) (r A , r B , in) • G« (r B , r A , in)] 

hu 2 f°° f°° 

— *\ / duu 4 a A (iu)a B (iu) e- lu/c / d<?<? 



x e 













p iV 
«JV J 



- 2 



/(0-ff(0 : 



6at/c 2 



X 2 

-»(0-p 



b N 



N 

bN 
k 2 

N 



J 2 (qX) ) (74) 



comes from the cross term of bulk and scattering parts 
[with f(x) and g(x) being defined by Eqs. (57) and (58), 
respectively, £ = c/ (lu), and Z = z B — z A ] , and 

U AB (r A ,r B ) = [ duu 4 a A (iu)a B (iu) 

27T Jo 

x Tr [G (1) (r A , r B , in) • G (1) (r B , r A , in)] 



x e 



64tt 3 7 

-(bw+b'jv)^ 



/>' r s r p ' 
U N N N 



duu a A (iu)oeB(iu) \ dqq I dq q 
o Jo 



r N r N 

b N b' N 



p jp 1 



r\ T r 



N' N 



k 4 
N 



bNb' N + 



2q 2 q 



2„f2 



bNr%r p N 



b N k- 2 



b' N k N ^ 



b N b' N 
J (qX)J (q'X) 



+ Aqq ':f N MqX)J 1 {q<X) + 



+ 



k 4 
N 

bNb' N r N r% 
k 4 



W r s r p ' 
N N N 



b N r%r 



N 1 N 

^Nb'jy 



JV 



b' M k* 2 



b N k N o n k n 

X J 2 {qX)J 2 (q'X) 



(75) 



is the scattering-part contribution [6^, = b^iq' , n), r' 
= r N (q',u)]. 



1. Perfectly reflecting plate 

Let us consider the case N = 1 (Fig. 2) in more detail 
and begin with the limiting case of a perfectly reflecting 
plate, 



r p = r p 



±1, 



(76) 



where the upper (lower) sign corresponds to a perfectly 
conducting (permeable) plate. In the retarded limit, 



where I, za, ZB»c/w min [u 
AO) 



min({^,|^' = A, J B; 



1,2,...})], U k Xb is g iven b Y E q. (62) with n(0) = 1 = 
n(0), whereas U^l [Eq. (74)] and U ( A 2) B [Eq. (75)] can be 
given in closed form only in some special cases. If X 
<C Z + (cf. Fig. 2), we derive, on using the relevant ele- 
ments of the scattering Green tensor as given in App. C 
[Eqs. (C10) and (Cll)], 



U 



23 l 3 Z+(l + Z+f 

(2) Cr 



AB 



(77) 
(78) 



where C r is given by Eq. (63) with e(0) = 1 = n(0). Thus, 
recalling Eq. (62), the interaction potential (72) reads 



U AB = 



1 32 X 2 + 6l 2 1 
f 23 PZ+(l + Z+f ~ ~Z\ 



a . (79) 
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A Z 



X 



App. C [Eqs. (C14)-(C17)], 



U 



(i) 



AB 



± 



4X 4 - 1Z 2 Z\ + X 2 (Z 2 + Z 2 ) 



C n 



U 



(2) 



AB 



(84) 
(85) 



= yX 2 + Z+), where C nr is given by Eq. (65) with 

e(m) = l. Hence, the interaction potential (72), reads, on 
recalling Eq. (64), 



FIG. 2: Two atoms near a perfectly reflecting plate. 



In particular, if za <C zb, or equivalcntly Z + ~ Z 
from Eqs. (77) and (78) it follows that 



(0) 

AB ' 



(0) 



u 



(2) 



(80) 
(81) 



so the interaction potential Uab, Eq. (72), is enhanced 
by the presence of the perfectly reflecting plate: 



U AB 



( 40^(0) 



52^(0) 

23 U AB 



for r„ fs1 = ( ±)1, 
forr p(s) =(+)!. 



(82) 



Next, we discuss the behavior of Uab in the case where 
the condition z^Czs is not valid. Since the bulk part 



U 



(o) 

AB 



[first term in the square brackets in Eq. (79)] is 



negative, the interaction potential is enhanced (reduced) 
by the plate if the scattering part U AB + U AB [second 
and third terms in the square brackets in Eq. (79)] is 
negative (positive). In the case of a perfectly conducting 
plate, it is seen that especially for Z = 0, briefly referred 



to as the parallel case, U AB + U AB is positive, and hence 
the interaction potential is reduced by the plate, whereas 
for X = 0, briefly referred to as the vertical case, U AB + 

(2) 

U\b i s positive and the interaction potential is reduced 
iff 



(2) 



zb/za < 4.90, 



(83) 



where, without loss of generality, atom A is assumed to 
be closer to the plate than atom B. It is apparent from 
Eq. (79) that for a perfectly permeable plate U AB + U AB 
is always negative, and hence the interaction potential is 
always enhanced by the plate. 

Let us now turn to the nonrctarded limit, where 



l,z A ,z B < c/w max [u 



max({w^,|A' = A, B; 



1,2,.. .})], and U A % is given by Eq. (64) [e(iu) = 1]. From 
Eqs. (74) and (75) we derive, on making use of the rele- 
vant elements of the scattering Green tensor as given in 



Uab = 



1 AX A 



X 2 (Z\ + Z 2 ) 1 



31% 



Cnr • 

(86) 

Let us again consider the effect of the plate on the 
interaction potential for the parallel and vertical cases. 
In the parallel case, Eq. (86) takes the form 



Uai 



1 



AV 



Z 2 



1 



I 6 3P(i 2 + z 2 + )i (i 2 + z 2 + r 

which in the on-surface limit Z + ^0 approaches 



C nr , (87) 



Uab 



2 f/ (0) 
3 U AB 



for r. 



P(s) - (-) 



f U AB imr p(s) 



(+) 



1. 



It can be seen easily that the term U AB [second term 
in the square brackets in Eq. (87)] dominates the term 
U AB [third term in the square brackets in Eq. (87)], so 

U AB + U AB is positive (negative) for a perfectly conduct- 
ing (permeable) plate, and hence the interaction poten- 
tial is reduced (enhanced) due to the presence of the 
plate. 

In the vertical case, from Eq. (86) the interaction po- 
tential is obtained to be 



Uab 



It is obvious that U 



1 2 

¥ T 3ZjP 



1 



(89) 



(i) 

AB ' 



r(2) 

AB 



Eq. (89)] is negative when the plate is perfectly conduct- 
ing, thereby enhancing the interaction potential since 

U AB [first term in Eq. (89)] is negative. In the case of a 

perfectly permeable plate, U AB + U AB is positive iff 



^<1 

ZA 



1 



14.82, 



(90) 



where atom A is again assumed to be closer to the plate 
than atom B. 



Since U AB and U AB arc negative in all cases, the re- 
alization of enhancement or reduction of the interaction 
potential depends only on the sign of U AB and its magni- 

(2) 

tude compared to that of U AB . In particular, the results 



r(2) 
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A' 



FIG. 3: Two electric dipoles near a perfectly conducting plate 
(parallel case). 



FIG. 5: Two magnetic dipoles near a perfectly conducting 
plate (parallel case). 





FIG. 4: Two electric dipoles near a perfectly conducting plate 
(vertical case). 



FIG. 6: Two magnetic dipoles near a perfectly conducting 
plate (vertical case). 



for the non-retarded limit (the sign of U^ B being sum- 
marized in Tab. I) can be explained by using the method 
of image charges, where the two-atom vdW interaction 
is regarded as being due to the interactions between fluc- 
tuating dipoles A and B and their images A' and B' in 
the plate, with 

Hint = Vab + Vab> + V BA > (91) 

being the corresponding interaction Hamiltonian. Here, 
Vab denotes the direct interaction between dipolc A and 
dipole B, while Vab 1 and Vba 1 denote the indirect inter- 
action between each dipole and the image induced by the 
other one in the plate. The leading contribution to the 
energy shift is of second order in H lnt , 

a Z7 W (°a,0b| .Hinting, m B ) 

AE AB --^ -j——— 

n,m \ si a / 

x (n A ,m B \H int \0 A ,0 B ). (92) 

In this approach, U A ° B corresponds to the product of 
two direct interactions, so it is negative in agreement 
with Eq. (86), because of the minus sign on the r.h.s. of 



Eq. (92). Accordingly, U AB is due to the product of two 
indirect interactions and is also negative — in agreement 
with Eq. (86). The terms containing one direct and one 

indirect interaction are contained in U^ B and determine 
its sign. We can hence predict the sign of U^ B from a 
graphical construction of the image charges, as sketched 
in Figs. 3-6. 

Figure 3 shows two electric dipoles in front of a per- 
fectly conducting plate in the parallel case. The configu- 
ration of dipoles and images indicates repulsion between 
dipole A(B) and dipole B'(A'), so U^ B is positive, in 
agreement with Tab. I. On the contrary, in the verti- 
cal case from Fig. 4 attraction is indicated, i.e., negative 

^AB' w hich is a l so in agreement with Tab. I. 

The case of two electric dipoles in front of a perfectly 
permeable plate can be treated by considering two mag- 





conducting plate 


permeable plate 


parallel case 






vertical case 




+ 



TABLE I: Sign of U^b f° r a perfectly reflecting plate. 
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nctic dipolcs in front of a perfectly conducting plate, as 
the two situations are equivalent due to the duality be- 
tween electric and magnetic fields in the absence of free 
charges or currents. From Figs. 5 (parallel case) and 6 
(vertical case) it is apparent that the interaction between 
dipole A(B) and dipole B'(A') is attractive in the parallel 
case and repulsive in the vertical case, again confirming 

the sign of U A X B as gi ven m Tab. I. When the dipole- 
dipole separation in Fig. 6 is sufficiently small compared 
with the dipole-surface separations, then the direct inter- 
action between the two dipoles is expected to be stronger 
than their indirect interaction via the image dipoles. As 
a result, U AB will be the dominant term in U AB + U AB 

and U AB + U^b becomes positive. However, when the 
dipole-dipole separation exceeds the dipole-surface sep- 
arations, then the indirect interaction may become com- 

(2) 

parable to the direct one, and U AB may be the dominant 



term, leading to negative U AB + U AB . The image dipole 
model hence gives also a qualitative explanation of the 
condition (90). 



r(2) 



2. Semi-infinite magnetodielectric half space 

Let us now abandon the assumption of perfect reflec- 
tivity and consider a magnetodielectric plate of permit- 
tivity e(u>) and permeability j-i(u>). To be more specific, 
we restrict our attention to a sufficiently thick plate so 
that the model of a semi-infinite half space applies. In 
this case, Eq. (71) for the reflection coefficients reduces 
to 



Kb 



bo 



(93) 



with b = bi = ^/u 2 /c 2 +q 2 , bo = y/s(iu)/j,(iu)u 2 /c 2 +q 2 , 
Ag=^(i?i), and Ag=e(m). 

In the retarded limit, /, za, zb c/cj m j n [with ui m i n be- 
ing defined as above Eq. (61)] we may again replace the 
atomic polarizability and the permittivity and permeabil- 
ity of the plate by their static values. Replacing the inte- 
gration variable q in Eq. (74) by v = b\c/u [cf. Eq. (C9)] 
leads to 



U^ta^b) 



he 



32tt s Ps 



IK 



r 

2«A(0)aij(0) / 
J\ 



Z 2 A 5 _ + (Z 2 - 2X 2 ) 



A A 



+ 



dv 

A 3 - 



+ l 2 A 5+ + IA 4+ + A 3+ + 2(v 



I 2 

1)( X 2 B 5 



+ (X 2 - 2Z 2 ) 



Bj B 3 
I + I 2 



Z 2 A 5+ + (Z 2 - 2X 2 



+ l 2 A 5 - + IAa- + A 3 . 



A 



4+ 



A 



3+ 



r- 



(94) 



where, according to Eq. (93), the static reflection coeffi- 
cients are given by 



.,,_MQ)t>-v^6y^H? (95) 



fi{o)v + v^oRo) - i + 



e { *) V -^ )m -l + V 2 
P PV e(0)v+^/e(0)n(Q)-l + v 2 

1 f°° 

A n ± = — l du u n e- au / c [J (Pu/c) ± J 2 (pu/ 
c Jo 



and 



B n 



r n+l 



duu n e- au / c J Q ((3u/c), 



(97) 
(98) 



with P = X\Jv 2 — 1 and a = I + vZ + (for explicit expres- 
sions of A n ± and B n , see App. D). Similarly, Eq. (75) 
reduces to 



hu 2 



U {2) = — 

64^ 3 c 2 



a A (0)a B {0) / dv dv' 



x | [r p r' p [3v 2 v /2 - 2(v 2 + v' 2 ) + 2] + r s r' s - r s r' p v' 2 
- r p r' s v 2 ^M + ^vv'\Jv 2 - l^v' 2 - \r p r' p M x 
+ (r s r' s + r p r' p v 2 v' 2 + r s r' p v' 2 + r p r' s v 2 )M^ (99) 
[r' a = r a {v')], where 
\J„ = / duu 6 e- (v+v ' )z + u/c J n (f3u/c)J n (/3'u/c) (100) 



(/?' = X\/v' 2 — 1), which can be evaluated analytically 
only in some special cases. In particular, when X <C Z+, 
then approximately 

M - J 2 <Q] [°° dvv 6 e-( v+v "> z + u / c - 720c7 S n 

M n - J n {U) ^ ciuu e - (y + d n0 . 

(101) 

In the nonretarded limit, I, za, zb -Cc/[n(0)a; max ] [with 
w max being defined as above Eq. (64)], U AB and U AB can 
be obtained by using in Eqs. (74) and (75), respectively, 
the relevant elements of the scattering part of Green ten- 
sor as given in App. C. In the case of a purely dielectric 
half space {fi = l) we derive [Eqs. (C20)-(C23)] 



Uai 



Cnr 



+ 



[4X 4 - 2Z 2 Z\ + X 2 {Z 2 + Z 2 )] C 



(i) 



(102) 



where C nr is given by Eq. (65) with e(iu) = 1, and 

h 



167r 3 eg J0 



3h 



e(iu) — 1 

duaA(iu)oeB(iu) , (103) 



l6^s 2 0J0 



du a>A(iu)a>B(iu) 



e(iu) + 1 ' 
e(iu) — 1 



1 2 



e(iu) + 1 



(104) 
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In particular in the limiting case when /<CZ+, Eq. (102) 
reduces to 



C nr ^ (X 2 -2Z 2 )c£ ) 

U AB - -— + ^3 



(105) 



It is seen that the second term on the r.h.s. of this equa- 
tion is positive (negative) in the parallel (vertical) case, 
so the vdW potential is reduced (enhanced) by the pres- 
ence of the dielectric half space. In the case of a purely 
magnetic half space (e= 1) we derive [Eqs. (C24)-(C27)] 



Uab = 
where 



Cnr 



+ 



647T 3 £ 2 |C 2 Ju 



[Z 2 - 2X 2 + 3Z+(l+ - Z+)]C, 



du u 2 aA (iu)a b (iu) 
fj,{iu) + 1 



(3) 



(106) 



(107) 



(2) 

Note that U AB does not contribute to the asymptotic 
nonretarded two-atom vdW potential Uab for the purely 
magnetic half space. In particular in the limiting case 
when X< Z+, Eq. (106) reduces to 



C7 nr (2Z 2 -X 2 )q 
Uab--—+ ^j-— 



(3) 



(108) 



It is seen that the second term in the r.h.s. of this equa- 
tion is negative (positive) in the parallel (vertical) case, 
so the vdW potential is enhanced (reduced) due to the 
presence of the magnetic half space. 

It should be pointed out that the nonretarded limit 
for the magnetodielectric half space is in general incom- 
patible with the limit of perfect reflectivity [s(iu] — > oo 
or jti(iti) — > oo] considered in Sec. IV Bl, as is clearly 
seen from the condition given above Eq. (102) [cf. also 
the expansions (C18) and (C19), which arc not well- 
behaved in the limit of perfect reflectivity]. As a con- 
sequence, Eq. (106) does not reduce to Eq. (86) via the 
limit fi(iu) — ► oo. It is therefore remarkable that the re- 
sult for a purely dielectric half space, Eq. (102), does re- 
duce to Eq. (86) in the limit e(iu) — > oo, as already noted 
in Ref. [39] in the case of the single-atom potential. 

Figures 7-9 show the results of exact (numerical) cal- 
culation of the vdW interaction between two identical 
two-level atoms near a semi-infinite half space, as given 
by Eq. (72) together with Eqs. (60), (74), and (75). In 
the figures the potentials and the forces are normalized 
w.r.t. their values in free space as given by Eq. (60) 
[n(iu) =1 = fj,(iu)]. In the calculations, we have used 
single-resonance Drude-Lorenz-type electric and mag- 
netic susceptibilities of the half space, 



e{u) = 1 



J Pe 



U), 



Te 



(109) 




FIG. 7: The vdW potential for two identical two-level atoms 
in the parallel case in presence of (a) a purely dielectric half 
space with wp e /a>io =3, WTe/wio = 1, and 7 e /u>io = 0.001 (b) 
a purely magnetic half space with /W10 = 3, WTm/oJW = 1, 
and 7 m /tuio = 0.001 is shown as a function of the atom- 
atom separation I [wio is the atomic transition frequency, and 
U AB (l) is the potential in free space]. The atom-half-space 
separations are za — zb = 0.01c/wio (solid line), 0.2c/o»io 
(dashed line), and c/wio (dotted line). 



MM = 1 + 



Tm 



(110) 



From the figures it is seen that the vdW interaction is 
unaffected by the presence of the half space for atom- 
half-space separations that are much greater than the in- 
teratomic separations, while an asymptotic enhancement 
or reduction of the interaction is observed in the opposite 
limit. 

Figure 7(a) shows the dependence of the normalized 
interaction potential Uab (I) on the atom-atom separa- 
tion I in the parallel case (Z — 0) for different values of 
the distance za (= zb) of the atoms from a purely dielec- 
tric half space. The ratio of the interatomic force along 
the connecting line of the two atoms, Fabx(1) [Eq. (55)] 
to the corresponding force in free space, F A °g x (l), follows 
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FIG. 8: The vdW potential for two two-level atoms in the ver- 
tical case in the presence of (a) a purely dielectric half space 
and (b) a purely magnetic half space is shown as a function 
of the atom-atom separation I. The distance between atom 
A (which is closer to the surface of the half space than atom 
B) and the surface is equal to O.OIc/cjio (solid line), 0.2c/wio 
(dashed line), and c/wio (dotted line). All other parameters 
are the same as in Fig. 7. 

closely the ratio Uab{1)/U^(1), so that, within the reso- 
lution of the figures, the curves for Fabx(Q / '^azL(0 ( not 
shown) would coincide with those for Uab(1)/U a %(1). 
The figure reveals that due to the presence of the di- 
electric half space the attractive interaction potential and 
force are reduced, in agreement with the predictions from 
the nonretarded limit, Eq. (105). The relative reduction 
of the potential and the force are not monotonic, there 
is a value of the atom-atom separation where the reduc- 
tion is strongest. The /-dependence of Uab(1) /U A g{l) m 
the presence of a purely magnetic half space in the par- 
allel case is shown in Figs. 7(b). Again, the correspond- 
ing force ratio Fabx(I) / F A g x (l) (not shown) behaves like 

UAB(l)/ufg(l). The figure indicates that the presence 
of a purely magnetic half space enhances the vdW in- 
teraction between the two atoms, with the enhancement 
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FIG. 9: The vdW force acting on atom A (which is closer to 
the surface of the half space than atom B) in the presence of 
(a) a purely dielectric half space and (b) a purely magnetic 
half space is shown as a function of the atom-atom separation 
I. All parameters are the same as in Fig. 8. 



increasing with the atom-atom separation, in agreement 
with the nonretarded limit, Eq. (108). 

Figure 8 shows Uab{1)/U a °1(1) in the vertical case 
(X = 0) when the half space is purely dielectric [Fig. 8(a)] 
or purely magnetic [Fig. 8(b)]. In the figure, atom A is 
assumed to be closer to the surface of the half space than 
atom B, and the graphs show the variation of the inter- 
action potential with the atom-atom separation I for dif- 
ferent distances za of atom A from the surface of the half 
space. It is seen that for a purely dielectric half space the 
potential is enhanced compared to the one observed in 
the free-space case — in agreement with Eq. (105). Note 
that there are values of the atom-atom separation at 
which the enhancement is strongest. For a purely mag- 
netic half space, the potential is seen to be typically en- 
hanced although for very small atom-atom separations a 
reduction appears [inset in Fig. 8(b)] — in agreement with 
Eq. (108). As in the parallel case, the relative enhance- 
ment is monotonic. Whereas the force Fbaz {1)/FjbAz (0 
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for the force acting on atom B (not shown) again fol- 
lows closely the potential ratio U AB (l) /U AB (l) > ^ ne ra tio 
Fab z (1)/F ABz (1), f° r tnc f° rce acting on atom A notice- 
ably differs from U AB {l)/U A a) B {l), as can be seen from 
comparing Figs. 8 and 9. Clearly, the reason must be 
seen in the different atom-atom and atom half-space di- 
rections in the two cases (cf. Figs. 4 and 6). 

Figs. 7(a) and 8(a) showing the interaction potential 
of two atoms in the presence of a purely dielectric half 
space in the parallel and vertical cases, respectively cover 
the results shown in Ref. [31] on a different scale. The re- 
sults here are more complete because they show that the 
relative potential does not have the monotonic behavior 
suggested by the figures in Ref. [31]. 

V. SUMMARY AND CONCLUSIONS 

Based on macroscopic QED in linear, causal media, we 
have obtained a general formula for the vdW potential of 
two ground-state atoms in the presence of an arbitrary 
arrangement of dispersing and absorbing magnetodielec- 
tric media by calculating the leading- (4th) order shift of 
the ground-state energy of the overall system. The result 
has been applied to two atoms (i) in bulk material (with- 
out taking into account local-field corrections), (ii) in the 
presence of a perfectly reflecting plate, and (iii) in the 
presence of a semi-infinite magnetodielectric half space. 
It has been found that the presence of a bulk magne- 
todielectric medium will reduce the interaction potential 
w.r.t. its well-known free-space value. 

We have further shown that in the presence of a per- 
fectly reflecting plate the vdW interaction can be en- 
hanced or reduced depending on the (electric/magnetic) 
nature of the plate and the (parallel/ vertical) alignment 
of the atoms. In particular, in the nonretarded limit these 
effects can be qualitatively explained using the method 
of image dipoles. 

Finally, we have calculated the vdW potential in the 
presence of a magnetodielectric half space. The analyt- 
ical results show that in the nonretarded limit the po- 
tential in the case of a purely dielectric half space is re- 



duced (enhanced) in the parallel (vertical) case compared 
to its value in free space, while in the case of a purely 
magnetic half space it is enhanced (reduced) for parallel 
(vertical) alignment of the two atoms. The results for a 
purely dielectric half space are in qualitative agreement 
with those for the perfectly conducting plate, while for 
a magnetic plate the results for finite permeability dis- 
agree with those for the perfectly reflecting case in the 
asymptotic power laws — owing to the fact that the two 
limits of perfect reflectivity and nonretarded distance do 
not commute. 

The numerical computation of the interaction poten- 
tial in the whole distance regime confirms the analyti- 
cal results. In addition, it shows that the relative en- 
hancement/reduction of the vdW interaction is not al- 
ways monotonous, but may in general display maxima or 
minima, in particular in the case of a purely dielectric 
half space. 
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APPENDIX A: INTERMEDIATE STATES AND 
INTERACTION MATRIX ELEMENTS 

The intermediate states contributing to the two-atom 
vdW interaction according to Eq. (34) are listed in the 
first three columns of Tab. II; the corresponding matrix 
elements of the interaction Hamiltonian (30) [together 
with Eqs. (8) and (25)] can be found by recalling the 
commutation relations (4) and (5) and using the relations 
(13) and (14). For example, for case (1) in Tab. II this 
leads to 



(1 (1) | (n A | (0s \H A f + H BF \0 A ) \0 B ) |{0}) 

= -[df-G* Xl (r A ,r 1 ,w 1 )]. i , (Al) 

(l {2 \l {3) \(0 A \(0 B \H AF + H BF \n A )\0 B )\l^) 
= -i={[d^.G^ ( r A ,r3,a;3)].^ 2 ) 

+ [dO".Gl 2 (r A ,r 2 , W2 )] j2 ^ 13 )}, (A2) 



(l^\(0 A \(m B \H AF + H BF \0 A )\0 B )\1 (2 \ 1< 3 )) 
= -i={[dr-G.3(rB,r3,c 3 )]^ (24) 

+ [dr-G A2 (r B ,r 2 ,^)] 42 «5( 34 )}, (A3) 

({0}\(0 A \(0 B \H AF + H BF \0 A )\m B )\l^) 

= -[d° B m -G X4 (r B ,r 4 ,w 4 )]. 4 , (A4) 
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Case |7) \II) \IH) Denominator 



(1) 


|n A ,0 s )|l (1) ) 


|0a,0 b }|1< 2 \i( 3 >) 


\Q A ,m B )\l^) 


A 


= (w3+w')(w' + w)(wb 










Ai 




(2) 


51 


\n A ,m B }\{0}) 


\0 A ,m B )\l^) 


Aii 


= M W)(^4+^b)(^B +W) 


(3) 


51 


51 


|nA,0 s >|l( 2 )> 


Av 


= (w2+u/)(w2 + wS)(w2+w) 


(4) 


n 


\n A ,m B )\l^,l^) 


\0 A ,m B )\l (4) ) 


A 


= (w3 + + w£ + u/ + a>)(w5 + w') 


(5) 


51 


51 


\n A ,0 B )\l w ) 


Am 


= (w5 + + ^ + ^' + w)(w3 + W) 


(6) 


\0 A ,m B )\l^) 


|0a,0 s )|1< 2 \i( 3 >> 


|nA,0 s >|l (4) > 


As 


= (w5+w')(w'+w)(wa+w'). 










Ai 


i = (wg + + + uj) 


(7) 


51 


\n A ,m B }\{0}) 


|nA,0 s >|l( 2 )> 


A x 




(8) 


51 


51 


\0 A ,m B )\l (2 >) 


A 


= (u# + w')(w5+wS)(u# + w) 


(9) 


51 


|n A ,m s >|l( 2 ),l( 3 )> 


|nA,0 s >|l (4) > 


Ai 


= (u# + w ')(w3 + + w' + + w') 


(10) 


51 


51 


\0 A ,m B )\l (4 >) 


Ai 


= (w£ + + + w' + + w) 



TABLE II: Intermediate states contributing to the two-atom vdW potential and corresponding denominators. 



where S^ a ^ is given by Eq. (37). Substituting them into 
Eq. (34), one obtains Eq. (36) and subsequently Eq. (39), 
with energy denominators A and Ai as given in Tab. II. 
The other denominators listed in the last column of the 
table follow in a similar way from the respective interme- 
diate states given in the first three columns. 



APPENDIX B: DERIVATION OF EQ. (40) 



From the energy denominators given in Tab. II, it is 
straightforward to obtain 



11 1 111 

- + ^ + — + ^ + ^ + A 



Ai Ai 
1 



1 



Ul^+U) Wr+W I \0) + w' uj — Ul' 



1 



Ul'X+W U'o + UJ' J \ U) + U)' UJ — UJ 



■ (Bl) 



Since the denominators appear in combinations of the 
form of Eq. (39), where they are multiplied with terms 
(the two factors in square brackets) which are always the 
same and symmetric with respect to w and uj', we may 
interchange uj <-> uj' in the second term and recombine it 
with the first one to obtain 



1 11 1 11 

7> + d- + d ••• + ~n ^ ~n ^ 77" 

J-^lll ^Vlll J-^IX -L^IV 



uj<X+uj% \ujI+uj 



1 



1 



UJi 



-UJ 



1 



UJ + UJ' 



UJ — UJ' 



(B2) 



where the symbol — > denotes equality under the double 
frequency integral. Similarly we have 



111 

A + 7a7 + 7a 



(u^+w')(wf \uj+uj' uj-uj' 

1 



(wg+w'XwX+wXw-a/) 



7w (B3) 



111 
n - 7 ~n~- n ■■ 

-Lsvn J^xi J^xn 

1 



(u>\+ uj'){uj™+uj') \uj+uy + uj-uj' 

1 

~ (w5 + u')(«b+w)(w-k/) 



(B4) 



The second terms in Eqs. (B3) and (B4) cancel each other 
after an interchange of u> <-» u/ to yield 



111 1 1 1 

a + £; + £~ + £~ + £~ + a^ 

2 / 1 1 



7 ■ (B5) 



(uJ^+U>)(u}g + Uj) \LU + Uj' UJ — UJ 



Summation of Eqs. (B2) and (B5) immediately leads to 
Eq. (40). 



APPENDIX C: SCATTERING GREEN TENSOR 
FOR THE PLANAR MULTILAYER SYSTEM 

The scattering Green tensor for a planar multilayer 
system can be given in the form [40] 

G (1) (r,r',m) = J d 2 qe lcl < r - r ">G (1) (q,z,z',iu) (CI) 
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(q _L e z ), where 



G (1) (q,z,z',m) = 



with 



8n 2 b N ^ e - e - r ^ e 

a— s,p 



sine 



COS( 



e p = Tt— (cos<pe x + sm0e y ) 



n 

k N 



(C2) 

(C3) 
(C4) 



(e 9 = cos <p e x + sin <f> e y = q/g, <z= |q|) denoting the polar- 
ization vectors for s- and p-polarized waves propagating 
in the positive(+)/ncgativc(— ) z-direction. Further, bn 
and fcjv, respectively, are defined according to Eqs. (69) 
and (70), and the generalized reflection coefficients are 
given in Eq. (71). Equations (C3) and (C4) imply that 



etel 



sin <p — sin <p cos < 

2 



sin (f> cos cos <j> 








(C5) 



cos 



fc-, 

ibivq 
k jv 



sm cos < 



COS( 



-#sind) cos < 

-^sin 2 ^) 

-i&2 s in<£ 



ibivq 

ibivq 
k N 

q 



COS l. 



sm ( 



(C6) 



Substituting these results into Eqs. (CI) and (C2), per- 
forming the ^-integrals by means of [41] 



f 

Jo 



dx e % 



' co&(yx) = 2iri v J v {x), 



and using the relation 

Ji(x) J (or) - J 2 (x) 



(C7) 



(C8) 



we arrive at the Eqs. (66)-(68). 

In the particular case of a perfectly reflecting plate in 
the retarded limit, it is convenient to replace the integra- 
tion variable q in Eqs. (66)-(68) in favour of v = b\c/u, 
i.e., q — \Jv 2 — lu/c [see Eq. (69)], and hence 



f 

Jo 



dv — ■ 

c 



(C9) 



For X <C Z + , the exponential terms effectively limits the 
integrals in Eqs. (66)-(68) to the region where qX <1, 
hence we can approximate J l/ (qX) by J„(0) = S^q, such 
that the nonzero scattering-Green tensor components 



read 

Gxx ( r a , r B , iu) = G$ (r A ,T B ,iu) 



1 



8irZ + 



1 



Z+u Z\_u< 



G^](r A ,r B ,iu) 
1 



2wZ + \ Z + u Z\u< 



r v e 



-Z+u/c 



-Z+u/c 



(CIO) 



(Cll) 



leading to Eqs. (77) and (78), recall Eq. (61). 

In the nonrctarded limit it can be shown that the main 
contribution to the frequency integrals comes from the 
region where u/{cb\) <C 1 (cf. Ref. [42]). In this region 
we have 



b\c 2 



h = b. 



(C12) 



By changing the integration variable q according to 

q 



J 

Jo 



dq 



bi 



/ db 

J u/c 



(C13) 



and setting the lower limit of integration to zero, from 
Eqs. (66)-(68) we find, after some algebra, the nonzero 
elements of the scattering Green tensor to be approxi- 
mately given by 



G^(r Al r B ,iu) 
G^j (r A ,r B , iu) 



2X 2 - Z\ 



Attu 2 



I 5 



1 pi 



Attu 2 l\ 



1 pt 



G i% x )^A,V B ,iu) = 

G { zJ ( r A,r B , iu) = 



_ c 2 3XZ+ 

(+) 47TU 2 1% 

2 X 2 -2Z\ 



Attu 2 1% 



1 pi 



(C14) 
(C15) 
(C16) 
(C17) 



with 1+ = ^X 2 + Z\, leading to Eqs. (84) and (85). 

For a semi-infinite magnctodiclcctric half space in the 
nonretarded limit, we apply a similar procedure as below 
Eq. (Cll) and expand the reflection coefficients given by 
Eq. (93) in terms of u/(bc), 

,2 



r s ~ 



fi(iu) — 1 /j,(iu)[s(iu)ii(iu) — 1] u 
H{iu) + 1 
e{iu) — 1 



\u{iu) + l] 2 b 2 c 2 
e{iu)[e{iu)u{iu) — 1] u 2 



e(iu) + 1 



[s(iu) + If 



(C18) 
(C19) 



Substituting (C18) and (C19) into Eqs. (66)-(68) and 
keeping only the leading-order terms of u/bc, in the case 
of the purely dielectric half space we can ignore r s and the 
second term in the r.h.s. of Eq. (C19), so the relevant el- 
ements of the scattering Green tensor are approximately 

m, x 2X 2 - Z 2 e(iu) - 1 c 2 

GV(r A ,r B , iu) = ^5 + eiu+1 ^. (C20) 
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G§}(r A ,r B ,iu) = 



1 e{iu) — 1 c 2 
Anil £ ( iu ) + 1 ^ ' 



(C21) 



APPENDIX D: EXPLICIT FORMS OF A n± AND 
B n IN EQS. (97) AND (98) 



Gxz( T A,rB,iu) = 



-G^J (r A ,r B , iu) 
3XZ+ e(iu) - 1 c 2 
" 4ttZ + e(m) + 1 u 2 ' 



(C22) 



. X 2 - 2Z 2 e(iu) - 1 c 2 

For a purely magnetic half space, the first term on the 
r.h.s. of Eq. (C19) vanishes, so the leading order of u/bc 
is due to the second term as well as the first term on 
the r.h.s. of Eq. (C18), so the nonzero elements of the 
scattering Green tensor can be approximated by 



l + — Z + fj,(iu) — 1 
AttX 2 [i(iu) + l 

, z + i + 



16ttX 2 ; + 



[H{iu) - 1], (C24) 



G$ (r A ,r B , iu) = 



+ 



16 ^M-i] 

Z+1+ — Z 2 ^ jj,{iu) 



4irX 2 l+ n(iu) + 1 ' 



(C25) 



G^](r A ,r B ,iu) = -G^J (r A , r B , iu) 
1+-Z+, 



16ttXI+ 



\p(iu) - 1], 



(C26) 



The integrals in Eqs. (97) and (98) can be performed 
to obtain the following explicit expressions: 



A. 



6a 



3+ 



14+ 



(fl 2 +/? 2 ) 5 

6 (a 3 -4a/? 2 ) 
(a 2 + /3 2 ) i 
6(4a 2 -/3 2 ) 



A 5+ = 



A 5 - = 



(a 2 +/3 2 ) 3 

6(4a 4 -27a 2 /3 2 + 4^ 4 ) 

(a 2 + /3 2 )^ 
30(4a 3 - 3a/? 2 ) 

(a 2 +/3 2 )^ 
30(4a 5 - 41a 3 /3 2 + 18a/3 4 ) 



B* = 



Ba = 



Bn = 



(a 2 + [3 2 ) 2 
3a(2a 2 - 3/3 2 ) 

(a 2 + /3 2 )' 
3(8a 4 -24a 2 /3 2 + 3/3 4 ) 

(a 2 +/3 2 ) f 
15a(8a 4 - 40a 2 /? 2 + 15^ 4 ) 
(a 2 +/3 2 )^ 



(Dl) 
(D2) 
(D3) 
(D4) 
(D5) 
(D6) 
(D7) 
(D8) 
(D9) 



G^J (r A ,r B ,iu) = 



16ttZ + 



[li(iu) - 1]. 



(C27) 
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